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ABSTRACT 



^3 



O . 

, Recent iV-body simulations show that the formation of a present-day, galaxy sized 

■ dark matter halo in the cold dark matter cosmogony in general consists of an early 

fast collapse phase, during which the potential associated with a halo is established, 
followed by a slow accretion phase, during which mass is added rather gently in the 
outer region. In this paper, we consider the implication of such a halo assemble history 
for galaxy formation. We outline a scenario in which the fast collapse phase is accom- 
panied with rapid formation of cold clouds and with starbursts that can eject a large 
amount of gas from the halo centre. The loss of orbital energy of the cold clouds to 
\ the dark matter and the ejection of gas from the halo centre by starbursts can signifi- 

^ ■ cantly reduce the halo concentration. The outflow from the starburst can also heat the 

^\ ' gas in the protogalaxy region. Subsequent formation of galaxies in the slow accretion 

l/^ \ regime is therefore in halos that have been pre-processed by these processes and may 

t;!" ' have properties different from that given by A''-body simulations. This scenario can 

help to solve several outstanding problems in the standard ACDM model of galaxy 
formation without compromising its success in allowing structure formation at high 
redshift. The predicted rotation curves can be significantly flatter than those based 
on the halo profiles obtained from iV-body simulations, alleviating the discrepancy of 
r"| . the TuUy-Fisher relation predicted in the standard ACDM model with observations. 

C^h' The flattened galaxy halos allow accreted minihalos to survive in their central regions 

longer, which may be helpful in producing the flux anomalies observed in some grav- 
^ ■ itational lensing systems. The preheating by the early starbursts effectively reduces 

^ ' the amount of gas that can be accreted into galaxy halos, which may explain why the 

^ , baryon fraction in a spiral galaxy is in general much lower than the universal baryon 

fraction, fs ^ 0.16, in the standard ACDM model. 
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1 INTRODUCTION 

For many years, the study of galaxy formation has to con- 
tend with uncertainties both in cosmology (which deter- 
mines the spacetime background and initial conditions for 
structure formation) and in the physical processes that 
drive evolution. Thanks to recent observations on the Cos- 
mic Microwave Background (CMB) anisotropy and large- 
scale structure, the cosmological uncertainties are largely 
resolved, and the study of galaxy formation can now focus 
on the physical processes. There is now much evidence that 
we live in a flat universe which is dominated by Cold Dark 
Matter (CDM), with a total energy density Q,o = 1.02±0.02, 
total matter density Hm.o = (0.135±0.009)/i^^, baryon den- 
sity f^B.o = (0.0224 ±0.0009)/l"^ Hubble constant (in units 
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of 100kms"^Mpc"^) h = 0.71 ± 0.04, the power-law in- 
dex of initial perturbation n ^ 1, and the amplitude of the 
perturbation power spectrum, as specifled by the rms of the 
perturbation field smoothed in spheres of a radius 8 /i~^Mpc, 
(Tg = 0.85 ±0.1. The values quoted here are based on the re- 
cent results obtained from WMAP observations of the CMB, 
but are consistent with a large number of other observations 
(e.g. Spergel et al. 2003) . Assuming the density perturbation 
to be Gaussian, this (small) set of parameters completely 
determines the CDM model for structure formation. Thus, 
although the nature of the dark matter and the cosmological 
constant remains to be explained, there is not much freedom 
in choosing the spacetime background and initial conditions 
for galaxy formation. 

When combined with a set of standard (simple) assump- 
tions about galaxy formation, this cosmogony has several 
outstanding problems. First, A-body simulations have con- 
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stantly shown that dark matter halos in this cosmogony are 
quite concentrated, which seems to be at odds with the ob- 
served rotation curves of dwarf galaxies (e.g. Moore 1994; 
Burkert 1995; van den Bosch et al. 2000; de Blok, McGaugh 
& Rubin 2001; Borriollo & Salucci 2001; Blais-Ouolcttc, Am- 
ram & Carignan 2001; Button et al. 2003) and the observed 
Tully-Fisher relation for normal disk galaxies (Mo & Mao 
2000; Navarro & Steinmetz 2000; van den Bosch et al. 2003). 
Because of this, alternative models are proposed, where ei- 
ther the dark matter is assumed to be warm or coUisional, or 
the initial power spectrum is assumed to be strongly tilted 
(e.g. Spergel & Steinhardt 2000; Hogan & Dalcanton 2000; 
Kamionkowski & Liddlc 2000). Although these models may 
help to reconcile some of the problems outlined above, they 
are not favored by the recent WMAP result that the format 
tion of the first generation of ionizing sources may have oc- 
cured at quite high redshift. The other problem is an old one. 
When modern CDM dominated theories of galaxy formation 
were first considered in the late seventies (e.g. White & Rees 
1978), it was immediately realized that some processes must 
be invoked to prevent gas from cooling too fast, otherwise 
all the gas would cool and turn into stars by the present 
day. Although various feedback schemes have been proposed 
(Dekel & Silk 1986) to solve this 'overcooling' problem, a sat- 
isfactory solution remains lacking. In the standard ACDM 
model favored by current observations, the problem is much 
more severe since the baryon fraction is approximately 16%, 
which is 3 times higher than that in the standard CDM 
model with f2m = 1, making radiative cooling even more ef- 
fective. The observed baryonic mass fraction in present-day 
galaxy halos is only a few percent (see §2.3); this poses an- 
other serious problem for the present cosmogony , because 
one has to explain why only such a small fraction of the gas 
in a protogalactic region eventually assembles into a galaxy. 
The conventional feedback scheme, where energy feedback 
from supernova explosions is assumed to drive the cooled 
gas in a formed galaxy disk out of a galactic halo, may not 
be a viable solution, because such feedback is only effective 
in very low mass galaxies (e.g. Mac Low & Ferrara 1999). 

The existence of these problems is not necessarily a fail- 
ure of the standard ACDM model because galaxy formation 
involves many complex physical processes in addition to the 
cosmological initial and boundary conditions. It is, therefore, 
possible that we are indeed living in a ACDM universe, but 
that some important pieces of the puzzle are still missing 
in our current understanding of galaxy formation. In par- 
ticular, the dynamical role of the baryon component and 
star formation in structure formation are two important as- 
pects of the problem that have not been fully understood. 
Since the universal baryon fraction in current cosmogony is 
quite high and since baryonic gas is dissipational, the baryon 
component may have played an important role in structure 
formation, especially on galactic scales. In this paper, we 
examine the potential of such dynamical effects in solving 
the problems discussed above. 

The scenario we are proposing is motivated by the fol- 
lowing considerations. According to recent A'^-body simula- 
tions (e.g. Wechsler et al. 2002; Zhao et al. 2003a,b), the 
formation of a galaxy halo in the ACDM cosmogony typi- 
cally consists of two distinct phases: a fast collapse phase 
during which the gravitational potential associated with a 
halo is established, and a slow-accretion phase during which 



mass is added gently to the halo without significantly af- 
fecting the potential well. The early fast-collapse phase is 
dominated by violent mergers of dark matter halos. Since 
radiative cooling of gas is very effective in galaxy-sized halos 
at high redshift (with cooling time much shorter than halo 
collapse time), this phase is expected to be accompanied by 
rapid formation of dense gas clouds, which then rapidly col- 
lapse into the centre of the galaxy halo. During this process, 
the structure of dark matter halos can be affected by the 
baryonic component. The violent collapse of the gas clouds 
in the fast-collapse phase may also trigger an episode of 
rapid star formation that may drive a large amount of gas 
out from the halo centre. Such mass loss will cause the halo 
to expand in the inner region, thereby reducing the halo 
concentration. As we will show below, these two processes 
combined may be sufficient to overcome the concentration 
problem in the standard ACDM model, given the high bary- 
onic mass fraction in this cosmogony. This early episode of 
starbursts associated with the fast-collapse phase may also 
heat the gas in a protogalaxy before it collapses, thereby 
reducing the amount of gas that can be accreted during the 
slow- accretion phase. A proper treatment of these processes 
is therefore also a key step in understanding the baryon frac- 
tion in normal galaxies (e.g. Mo & Mao 2002; Oh & Benson 
2003; Granato et al 2003). 



Some of these issues we are considering here have al- 
ready been discussed in earlier investigations. The impor- 
tance of the baryon component in affecting halo profile was 
pointed out in Binney et al. (2001); El-Zant et al. (2001) 
considered the possibility that the interaction between gas 
clouds and dark halo can erase the cusps in the halo of dwarf 
galaxies; Weinberg & Katz (2002) found that the secular 
evolution of a bar-like structure in a dark matter halo can 
effectively get rid of the cusps of CDM halos; Navarro et 
al. (1996) considered the possibility that a sudden blowout 
of baryons from the centre of a dwarf galaxy can cause sig- 
nificant flattening of the inner profile of its halo. The effect 
of preheating on subsequent gas cooling and galaxy forma- 
tion has been considered in Mo & Mao (2002) and in Oh 
& Benson (2003). Our study here intends to present a re- 
lated discussion in the context of dark halo formation in 
the current CDM cosmogony, and in particular, we want to 
demonstrate that a proper treatment of these processes may 
be pivotal in solving many of the vexing problems in current 
models of galaxy formation. 



This paper is organized as follows. In Section 2 we re- 
visit some of the problems about galaxy formation in the 

standard ACDM model. In Section 3 we examine two possi- 
ble ways in which the baryon component can affect halo con- 
centration, and we discuss the observational consequences of 
the change of halo concentrations in Section 4. In Section 5, 
we examine how an early phase of starburst, which may be 
responsible for the formation of bulge and halo stars, can 
affect the subsequent formation of galaxies in dark halos. 
Finally, in Section 6, we make further discussion about our 
results. 
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2 SOME OUTSTANDING PROBLEMS OF 

GALAXY FORMATION IN THE STANDARD 
ACDM MODEL 

2.1 The properties of CDM halos 

One important aspect of gala^xy formation in the CDM cos- 
mogony is the formation of dark matter halos. Since dark 
matter particles take part only in gravitational interaction, 
TV-body simulations can be used to study in detail the for- 
mation of dark halos in the cosmic density field, and a great 
deal has been learned about the properties of the halo pop- 
ulation in the past few years. One important finding is that 
CDM halos can be approximated by a universal profile of 
the form 

4ps 



p(r) 



(1) 



(r/rs)(l-Fr/rs)2 ' 

where rs is a scale radius and ps is the density at this rar 
dius (Navarro et al. 1997, hereafter NFW; see Moore et al. 
1998, Jing 2000, Klypin et al. 2001 for discussions about the 
uncertainty of the inner profile). The proper size of a halo 
is often defined so that the mean density within the halo 
radius rh is a factor Ah times the mean density of the uni- 
verse, p, at the redshift z in consideration. The halo mass 
Mh is then related to the halo radius by 



Mh 



-pAhrh . 



(2) 



In our discussion we adopt Ah 



(IStt^ -I- 82a; 



39a: )/i}in{z), where x — Qin{z) — 1, and ^1x1(2) is the mass 
density parameter at redshift z (Bryan & Norman 1998). 
For the standard ACDM model. Ah » 340 at z = 0. The 
parameter that characterizes the concentration of mass dis- 
tribution in a halo is the concentration parameter, defined 
to be c = rh/rs. It is then easy to show that the total mass 
within a radius r is 



Mh(< r) = M^f{cx)/f{c), 
where x = r/r\i, and 
f{x) = \n{l + x)-x/{l+x). 
The mass within rs is 
ln2-l/2,^ 
/(c) 



(3) 
(4) 

(5) 



Defining the circular velocity of a halo as Vk = (GMh/rh)^^^, 
we have 



Mh 



[Ahnm(z)/2]i/2Giy(z) 



(6) 



where H{z) is the Hubble constant at z. The gravitational 
potential of the halo at a radius r from the centre is 



1 ln(l -h cx) 



(7) 



/(c) X ' 

and the escaping velocity is Vescir) = y'— 2$(r). 

With large iV-body simulations, it is found that the 

typical value of the concentration depends on halo mass. 
For the standard ACDM model, this mass dependence can 
be parameterized as 

0.15 



c(M) = 11 



Mh 



1012/1-1 Mp 



(8) 



for halos identified at z = (e.g. Bullock et al. 2001; Zhao 
et al. 2003b). The value of c also depends on redshift. Nu- 
merical simulations show that the concentration of a halo is 
tightly correlated with its mass accretion history (Wechsler 
et al. 2002; Zhao et al. 2003a,b). As shown in Zhao et al. 
(2003a), the formation of each CDM halo in general con- 
sists of two phases: a fast collapse phase during which the 
gravitational potential is established (the formation of the 
core of potential) and a slow-accretion phase during which 
mass is gently added to the core of potential without sig- 
nificantly affecting the potential well. The formation of the 
core of potential is dominated by rapid mergers and halos in 
this phase have relatively low concentration, c ~ 4. In the 
slow- accretion phase, the values of inner quantities, such as 
Ts, Ms and Vs, change little, while c increases as more mass 
is added in the outer part of the halo. Based on this, one 
may define a formation time at which the circular velocity 
as a function of time remains roughly constant. Zhao et al. 
found that the formation time defined in this way is closely 
related to the concentration parameter, and for present-day 
halos the correlation is well described by 



'.Ax 



H{z() 



1 i/v 



Ha 



(9) 



where 77 ~ 1. We can invert this relation to write the forma- 
tion time as a function of c. For a flat universe, we have 



l + Zi: 



(3) 



2r, 



1/3 



-1/3 



(10) 



For halos of a given mass (or a given circular velocity), the 
median concentration defines a typical formation time tf, 
while the scatter in c implies a scatter in tf. The median 
of c for a \Q^^h~^ M© halo at the present time is about 11, 
corresponding to a formation redshift Z{ = 2. For a given 
mass, the dispersion in c is about 40%, which corresponds 
to a formation redshift between 1 and 3 for present-day halos 
with Mh ~ IQi^/i-^Mo. 



2.2 The concentration problem 

It has been pointed out by several authors that the cuspy 

NFW halos are at odds with the observed rotation curves of 
low-surface brightness galaxies, which are fit better by pro- 
files with a constant density core (e.g. de Blok et al. 2001). 
Even if one enforces the NFW profile in the fit, the obtained 
concentration parameters are usually much lower than that 
predicted by the standard ACDM model (e.g. Moore 1994; 
Burkert 1995; Borriello & Salucci 2001; Blais Ouelette, Am- 
ram & Carignan 2001). Although there have been concerns 
over the effects of beam smearing (van den Bosch et al. 2000; 
Swaters, Madore & Trewhella 2000), recent analyses of high- 
resolution Hq-HI rotation curves (dc Blok, McGaugh & Ru- 
bin 2001; and McGaugh, Barker & de Blok 2003) suggest 
that the discrepancy appears to be real. Similar problems 
may exist also for normal spiral galaxies. Matching model 
predictions with the observed TuUy-Fisher relation. Mo & 
Mao (2000) find that c ~ 4 is required to obtain an agree- 
ment. The constraint here comes from the fact that a higher 
halo concentration leads to a larger boost of the disk maxi- 
mum rotation velocity relative to the halo circular velocity, 
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and too high a concentration leads to too low a luminos- 
ity for a given maximum rotation velocity. In what follows, 
wo rc-visit this problem with the use of more recent data, 
putting it in the context of our discussions in the present 
paper. 



2.2.1 Problem with the Tully-Fisher Relation 

Spiral galaxies arc observed to show a tight relation between 
luminosity and rotation velocity Vobs (which is usually taken 
to be the maximum rotation velocity). This Tully-Fisher 
(TF) relation can be written in the form 



L = AV^, 



(11) 



where A is the TF amplitude and a the TF slope. For a 
given rotation velocity, the observed scatter in L is quite 
small, typically about 20% at the bright end. 

Theoretically, the mass of a halo with circular velocity 
14 is given by equation (6). If we assume the disk mass is 
Md = TTidMh, we can write 



Md = 9.3 X 10'" /i" 



M, 



200 kms-i 



/ rwd \ 
^ V 0.057 



200 



-1/2 



/v 



(12) 



where /v = Vobs /Hi is the boost of the observed rotation 
velocity relative to the circular velocity of the halo. We cal- 
culate the boost factor as a function of halo concentration c 
and disk mass fraction md, assuming that the halo response 
to disk formation is through adiabatic contraction (see e.g. 
Mo, Mao & White 1998 for detailed descriptions). The disk 
mass can be converted into a disk luminosity by assuming 
a disk mass-to-light ratio Td, and so equation (12) can be 
compared with the observed TF relation, once models for Td 
and /v arc adopted. Since /v depends on c, the observed TF 
amplitude is a constraint on c. This is how Mo & Mao (2000) 
obtained c ~ 4. 

One uncertainty in the predicted TF relation is the 
adopted mass-to-light ratio, which is assumed to be a con- 
stant for a given Vobs in the model while in reality it may 
vary from galaxy to galEixy. Thus, a better way to constrain 
the theoretical model is to use the observed baryon TF re- 
lation, which directly relates the mass of cold gas in a spi- 
ral galaxy with its rotation velocity Vohs- Such relation has 
been obtained by McGaugh et al. (2000) and Bell & de Jong 
(2002). By modeling in detail the stellar mass and cold gas 
mass in each galaxy in their sample, Bell & de Jong obtained 
a baryon mass of Md ?a 7.0 X 10^° Mo (assuming h = 0.71) 
for galaxies with Vobs = 200 km s^^, which is similar to the 
value, Md « 6.4 x 10^° M©, obtained by McGaugh et al. The 
typical scatter in Md at the massive end is about 20%. These 
observed amplitudes can be compared directly with the pre- 
diction given by equation (12), once we know the boost fac- 
tor /v (see e.g. Mo, Mao & White 1998). We also assume 
disks to have exponential mass profiles and the specific an- 
gular momentum is given by a spin parameter A = 0.04. As 
shown in Mo et al. (1998), this is the typical spin parameter 
required to reproduce the observed disk sizes. 

Fig. 1 shows the predicted TF amplitude at Vobs = 
200kms~^ as a function of rud for several values of c. A 
comparison with the observational results clearly shows that 
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Figure 1. (a) The predicted amplitude of the baryon TF rela- 
tion as a function of md and halo concentration c (assuming NEW 
profile) for disk galaxies with Vjbs = 200 km s^-'-. The solid hori- 
zontal line is the mean value observed by Bell & dc Jong (2001), 
while the two dashed horizontal lines represent 20% deviation 
from the mean, (b) The predicted amplitude of the B-band TF 
relation for V)bs = 200kms~^ as required to match the observed 
luminosity function. The horizontal solid line is the observed am- 
plitude based on TuUy & Pearce (2000), while the two dashed 
horizontal lines represent 20% deviation from the mean. 



a low concentration, c 5, is preferred. For c = 12, the me- 
dian value for CDM halos in the standard ACDM model, 
the predicted amplitude is lower by a factor of about two. 

Note that even if a low value of c is adopted, consistency 
with observation requires md to be substantially smaller 
than the universal value 0.16. Thus, only a relatively small 
fraction of baryons in a protogalgixy can manage to settle 
into the final disk. We will come back to the implication of 
this result in subsection 2.3. 

Recently Eke, Navarro & Steimetz (2001) revisited the 
zero-point of the TF relation using high-resolution numeri- 
cal simulations. They found that their zero-point is fainter 
than the observed value by about 0.5 magnitude for a 
fixed circular speed. They further argued that this dis- 
crepancy is not significant as their simulated galaxies have 
colours that are slightly too red compared with the ob- 
serve TF galaxies. Their result is therefore inconsistent with 
what we find here. In reaching their conclusions, however, 
they used circular velocities measured at a radius given by 
20(Vh/220kms^^)/i^^ kpc, instead of the maximum circular 
velocity usually adopted in observations. To show the effect 
of their particular choice of the circular velocity, we show in 
Fig. 2 the TF relation for md = 0.16 and c = 12; these values 
are similar to those in Eke et al. (2001). We also adopt an 
/-band mass-to-light ratio of l.lh. Results are shown for two 
A values (0.02 and 0.04). Their simulations likely correspond 
to the low A value as their galaxies have quite compact sizes. 
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Figure 2. The dashed and thin soUd lines are the predicted TuIIy- 
Fisher relation for A = 0.02 and 0.04, respectively. For each spin 
parameter A, the upper curve shows the TuUy-Fisher relation 
where the velocity corresponds to the maximum circular veloc- 
ity of the disk, while in the lower curve the circular velocity is 
measured at 20(yh/220 km s-l)/i-l kpc (as in Eke et al. 2001). 
The thick solid line shows the observed TuUy-Fisher relation by 
Giovanelli et al. (1997) in the 7-band. 



With the circular velocity definition adopted by Eke et al. 
(2001), the predicted TF relation matches well the observed 
TF relation. This confirms their result. However, if we use 
the maximum circular velocity, as in observations, the pre- 
dicted TF amplitude is fainter than the observed one by 
3.0 and 1.0 mag, respectively. This exercise again highlights 
that with a large and a large concentration parameter, 
the TF zero-point is difficult to match. 



2.2.2 Problem with the Galaxy Luminosity Function 

Another constraint on halo concentration can be obtained 
from the observed luminosity function of galaxies. Detailed 
modeling of the luminosity function of galaxies in the 
standard ACDM by Yang et al. (2003; see also van den 
Bosch et al. 2003b) shows that the halo mass-to-light ra- 
tio is about lOO/i (in the _B-band) for halos with masses 
Mh ~ lQ^^h~^ M0 and is higher for both larger and smaller 
masses. This constraint comes from the fact that the number 
density of galaxy-sized halos is fixed in the standard ACDM 
cosmology, and so the observed luminosity density in the 
universe is directly related to the halo mass-to-light ratio. 
Note that this result of the mass-to-light ratio for galaxy 
halos is also consistent with the result obtained with grav- 
itational lensing (e.g. Hoekstra, Yee & Gladders 2003). For 
a given halo mass-to-light ratio, we can use equation (6) to 



write the the luminosity of the galaxy hosted by the halo as 

V 200 ) V200kms-iJ ' ^^'^> 

This has a form similar to the TF relation, with amplitude 
depending on the boost factor /v . The observed B-band TF 
amplitude at Vohs = 200 kms"^ is 8.6 x lO^/i"^L0. This 
value is quoted without including dust extinction correc- 
tion, because such correction is not made in the mass-to- 
light ratio quoted above (see Yang et al. 2003 for detailed 
discussion). The constraint given by this relation is shown 
in the lower panel of Fig. 1. To predict a high enough lu- 
minosity with c = 11, the m,^ is required to be extremely 
low, but such a low is not allowed by the baryon TF rela- 
tion. This inconsistency between the TuUy-Fisher zero-point 
and the luminosity function in the standard ACDM model 
has also been found in almost all semi-analytical models of 
galaxy formation (e.g. Sommerville & Primack 1999; Kauff- 
mann et al. 1999; Benson et al. 2002; Mathis et al. 2002; van 
den Bosch, Mo & Yang 2003a). As one can see from the fig- 
ure, consistent results can be obtained when c is lower than 
the typical CDM value. 



2.3 Baryon fraction in spiral galaxies 

Figure 1 shows that consistency with observation requires 
rrid to be much smaller than the universal value 0.16, even 
if a low value of c is adopted. Thus, only a small fraction of 
baryons in a protogalaxy can manage to settle into the final 
disk (see also the discussion in Salucci & Persic 1997). A sim- 
ilar conclusion can be reached from detailed modeling of the 
mass components in the Milky Way. Based on the kinemat- 
ics of Galactic satellites and halo objects, the halo mass of 
the Milky Way is estimated to be about 2 x 10^^ M0 within a 
radius of about 200 kpc (Wilkinson & Evans 1999; Sakamoto 
et al. 2002). Within this radius, the mass is still dominated 
by the Milky Way halo, because it is much smaller than 
the distance to the other bright galaxy in the Local Group, 
M31. The K-band absolute magnitude of the Milky Way ob- 
tained by Drimmel & Spergel (2001) is about —24.02, cor- 
responding to 8.3 X IO^^Lq. Assuming {M/L)k ~ 1 (e.g., 
Binney & Merrifield 1998), we find a total stellar mass of 
about 8.3 X 10^" Mq, which gives a baryon/total mass ratio 
of ~ 4.2%. A similar number can be obtained from the local 
baryonic surface density, which is approximately 50 M0 /pc^ 
(Kuijken & Gilmore 1991). Assuming the Galactic Center is 
8 kpc away, and taking a disk scale-length of 2.24 kpc (Drim- 
mel & Spergel 2001), one obtains a mass 5.6 x 10^° Mq for 
the disk. The bulge/disk mass ratio is about 1/5 (Kent et 
al. 1991), and so the total mass in the bulge and disk is 
about 6.7 X 10^" Mq, which gives a baryon/total mass ratio 
of ~ 3.4%. Since the mass of baryons in other components, 
such as cold and hot gas, is much smaller, the total baryon 
fraction in the Milky Way is therefore much smaller than 
the universal value 16%. 
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3 DYNAMICAL EFFECTS OF THE BARYON 
COMPONENT 

The problems presented above are not necessarily the failure 
of the standard ACDM model because, as mentioned earlier, 
galaxy formation involves many complex physical processes. 
One simple assumption made in the model predictions con- 
sidered above is that the growth of a galaxy disk in the halo 
centre is a gentle process, so that the gravitational effect of 
the baryon component can be modeled by adiabatic contrac- 
tion. This may be a good assumption, at least for halos with 
cuspy inner profiles, as shown by numerical simulations (e.g. 
Jesseit, Naab, Burkert 2002). In addition, the model also 
makes an implicit assumption that the density profiles of 
dark matter halos, within which disks grow, have the same 
properties as the halos in cosmological A^^-body simulations. 
This may not bo a good assumption if the formation of a 
galactic halo contains a fast-collapse phase during which the 
interaction between the baryon component and dark matter 
can significantly affect the structure of the halo. In this case, 
the formation of a galaxy disk is in a pre-processed halo that 
may have properties different from a CDM halo. 

Based on this consideration, we propose a scenario of 
galaxy formation that consists of the following aspects: 

(i) The formation of the core of potential through a fast 

accretion phase, during which gas cools and collapses to form 
dense clouds that interact dynamically with dark matter 
particles; 

(ii) The loss of orbital energy of dense clouds due to dy- 
namical friction against dark matter particles, which in- 
creases the kinetic energy of the dark matter in the inner 
region; 

(iii) A phase of rapid star formation in the fast-accretion 
phase, which causes the loss of a large amount of gas from 
the halo centre, and causes the inner part of the halo to 

expand; 

(iv) The heating of the gas in the protogala^xy region by 
the starburst in the fast accretion phase; 

(v) Subsequent formation of a galaxy disk in the slow 
accretion phase by cooling of the preheated gas in the pre- 
processed halo. 

In this section we focus on the effect of the baryonic 
component on the halo. 



3.1 Collapse of the baryon component and its 
interaction with dark matter 

In the fast collapse phase, i.e. during the formation of the 

core of potential, the interaction between the gas compo- 
nent and dark matter is a complicated process. In order to 
gain some insight into this process, let us first examine some 
relevant time scales. 

According to Zhao et al. (2003a), the mass of a halo in 
the fast collapse phase increases roughly as Mi^iz) oc (1 + 
z)~'^ in the standard ACDM model. For galax;y-sized halos 
where the fast collapse phases ended at this relation 

defines an accretion time scale, 

tac = M/M ~ l/[^H{z)] , (14) 
where we have used the fact that the scale factor increases 
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Figure 3. The solid curve shows the circular velocity as a func- 
tion of redshift for halos where cooling time is equal to the collapse 
time, while the dashed curve shows the masses of such halos. The 
cooling time uses the cooling function given by Sutherland & Do- 
pita (1993) for an ionized gas with metallicity equal to 0.1 times 
the solaj: value. 

with time as a{t) oc t^^^ at « ^ 2. This time scale is slightly 
shorter than the collapse time scale of the halo, 

tcoii ~ ^ ~ 1/[3H(^)] . (15) 

Note that if 14 is fixed, then - 2/[ZH{z)\. The fact that 
the accretion time scale in the fast collapse regime is much 

shorter than this is because the potential well associated 
with the halo deepens rapidly in this regime, with Vh oc 
{l + z)-^'\ 

If radiative cooling is effective, then gas will cool to form 
clouds before it collapses in the halo centre to form stars. 
For a completely ionized gas, the cooling time scale can be 

estimated as 



^cool 



(16) 



2n2A(T) ' 

where T is the virial temperature of the halo, n^K{T) is the 
cooling rate per unit volume, and n is the mean particle 
density. We use the cooling function given by Sutherland & 
Dopita (1993) for an ionized gas with metallicity equal to 
0.1 times the solar value. This should be compared to the 
collapse time fcoii • For halos with tcooi "C tcoii , gas cannot be 
heated up by the gravitational collapse, while for halos with 
icooi S> tcoii, gas may be heated to the virial temperature 
and form a hot halo which is approximately in hydrostatic 
equilibrium in the dark halo potential well. Fig. 3 shows the 
critical halo circular velocity and mass as a function of red- 
shift, calculated by equating tcooi and tcoii- This mass is quite 
independent of redshift and has a value about lQ^^h~^ M©. 
Thus, in the absence of heating, radiative cooling is effective 
for halos with masses below that of the Milky Way halo at 
all rcdsliifts. 

Because of the existence of perturbations on small 
scales, the collapse is expected to be very clumpy in the fast 
collapse phase of the halo. This, together with the rapid gas 
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cooling discussed above, implies that the gas associated with 
the collapse of galactic halos must be in the form of mas- 
sive cold clouds. Because of effective cooling, the gas clouds 
can acquire a high binding energy, so as to survive the tidal 
disruption by the dark halo. Once such gas clouds become 
self-gravitating, star formation may occur in situ before the 
coUpase of the final halo. However, based on the results of 
Kennicutt (1998), only about 1.7% of the gas in a galaxy (be 
it a normal spiral or a starburst) is converted into stars dur- 
ing a typical rotation time of the galaxy. Thus, the fraction 
of the gas that can be converted into stars during a collapse 
time of the halo is about 0.017(7rrh)/(27rrg), where rg is the 
radius of the galaxy. Since the gas typically contracts by a 
factor of 20-30 before it is supported by angular momen- 
tum, only a small fraction of the gas can be converted into 
stars during the fast collapse phase. Thus, much of the cold 
gas can collapse towards the halo centre. This is consistent 
with the fact that a lot of cold gas seems to be able to sink 
to the centre of a starburst galaxy to feed the intensive star 
formation there. 

As the cold clouds sink towards the centre of a dark 
matter halo, they interact with the dark matter through 
gravitation. We can obtain some ideas about the gravita- 
tional interaction between the cold clouds and dark matter 
particles from numerical simulations of major mergers of 
galaxies (e.g. Hcrnquist 1993). In such simulations, each of 
the merger progenitors is assumed to be a galaxy system 
consisting of a disk (sometimes a bulge is also iuchidcd) cm- 
bedded in an extended dark matter halo. Numerical simula- 
tions show that, when two such systems merge, their halos 
merge first to form a common halo, while the two galax- 
ies lose their orbital energy due to dynamical friction, and 
sink towards the halo centre to merge. Here the dynamical 
friction plays a key role in the merger of the galaxies. The 
situation in our problem is more complicated, because the 
baryon component is still in gaseous form, and so its energy 
can bo dissipated through cloud-cloud collisions. However, 
simulations of major mergers containing gaseous disks (e.g. 
Barnes 1992) suggest that such dissipation is important only 
when the gas clouds have already sunk to the centre of the 
halo, while most of the orbital energy of the cold gas may 
be lost through dynamical friction. In what follows, we use 
simple arguments to gain a qualitative understanding of the 
problem. 

According to Chandrasekhar's dynamical friction for- 
mula, the dynamical fiction force on a gas cloud with mass 
M is 



„ 47rG^(lnA)peffM% ^ 

r M = TTs » M ; 



(17) 



where Vm is the velocity of the gas cloud relative to the 
dark matter, InA ~ 10 is the Coulomb logarithm, and pcff 
is the effective mass density of halo particles producing the 
dynamical friction. The rate at which the gas cloud loses its 
orbital energy due to dynamical friction is 



E = Vm • Fn 



(18) 



For a gas cloud moving in a halo with density profile p{r), 
the energy loss rate at radius r is approximately 



For circular orbits, the decay of the orbital radius is given 
by 

dirVM) 



dt 



: rFu/M (X rpM/Vii . 



(20) 



Assuming that the energy loss of the gas cloud at radius r 
is deposited in a spherical shell with radius r ^ , the change 
in the specific energy of dark matter particles at radius r is 
then 



£ (X 



E 



4nr^p{r) dr 



dt M 

cx — y(r) , 



(21) 



where 



p{< r) 
p{r) 



1 + 



dlnFh 



dlnr 



p(< r) is the mean density of dark matter within radius r, 
Vii(r) = [GMh(< r)/rY^'^, and the last relation uses eqs. 
(19) and (20). 

The total energy loss of a cloud is 



EoTh ■ 



Jo 



4:npr^£{r)dr . 



(22) 



For an NFW profile, most of the energy loss is in the inner 

region, close to r — ts, where p oc r^^. The energy gain per 
unit mass by the dark matter at radius r is therefore 



Mo r 

where forb = Eoib/M, and 



Mo=g 



nr 



A'Krp{r)Q{r)dr . 



(23) 



(24) 



Over a large range of r. Mo is a constant, and Mo ~ Mh for 
an NFW profile with c = 4. 

As one can see, the effect of dynamical friction de- 
pends on the mass of gas cloud. If we approximate dark 
halo by a singular isothermal sphere, the dynamical friction 
time scale is approximately t-av ^ 0.3/(e)tcoii'?/ lii(')), where 
77 = M\i/M is the ratio between the halo mass and the cloud 
mass, /(e) ~ t^'"^^ describes the dependence on orbit eccen- 
tricity (Lacey & Cole 1993), and the typical value for /(e) 
is about 0.5 according to Tormcn (1997). Thus, clouds with 
masses M Mh/20 can lose most of their orbital energy 
through dynamical friction. During the fast collapse, the 
merger progenitors all have comparable mass (major merg- 
ers), and so dynamical friction is expected to be important 
for most cold clouds. 

Assuming isotropic velocity dispersion in a halo, the 
phase space density profile can be written as 



Q{r) 



pjr) 



(25) 



E<xM^p{r)/VMir). 



(19) 



where ai{r) is the original velocity dispersion at radius r. As 
shown by Taylor & Navarro (2001), this phase space density 
profile can be well approximated by a power law Q{r) oc 
^-1.825 Jqj. ]S[F"W density profile. If we assume the energy 
transfer from gas clouds at a radius r is to be thermalized 
completely around this radius, the effect on dark matter 



t This assumption holds only as an average for clouds moving on 
different orbits. 
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Figure 4. The phase-space density profile of dark matter parti- 
cles. The thick solid curve show that for the NFW profile, while 
the other solid curves show the profiles taking into account the 
effect of dynamical friction. Curves from bottom up arc results 
where the gas mass is 0.16, 0.08, 0.04 and 0.02 times the halo 
mass. The dashed curves show the results for the density profile 
in eq. (27) assuming c' = rj^/rc = 2.2, 3.5, 5.5, 8.0 (from bottom 
up). The values of Q at the outer radius of the halo are chosen 
to match those given by the solid curves. The vertical line marks 
the value of Ms (cf . eq 5) . 



particles is to increase their velocity dispersion and reduce 
the phase space density: 

p(r) 



Q{r) 



(26) 



[o-2(r)-FAo-2(r)]3/2 ' 

where A(T(r) is the increase in the velocity dispersion due 
to dynamical friction. Using the expression of £ to obtain 
Acr^, we can obtain the phase-space density profile. Figure 
4 show the results for different values of /gas = Mgas/Mh, 
assuming c = 4. As one can see, the phase-space density 
in the inner part of the halo can be significantly reduced, 
and the new phase-space density is roughly a constant in 
the central region. The effect is larger for a larger value of 
/gas- If /gas is comparable to the universal value, 0.16, more 
than 20% of the dark matter particles will have a roughly 
constant Q. 

There are a number of uncertainties in the results shown 
above. First of all, although Chandrasekhax's formula was 
shown to describe reasonably well the decay of a satellite's 
orbit in a halo of dark matter particles (e.g. Bontekoe & 
van Albada 1987; Zaritsky & White 1988; Weinberg 1986), 
the assumption of a local response of the halo may not be 
vahd (e.g. Weinberg 1989; Weinberg & Katz 2002). Because 
of this, our assumption that the decay of orbital energy of 
the gas clouds is deposited locally is clearly invalid in detail. 
Unfortunately, as discussed in Weinberg & Katz (2002), a 
reliable result of the halo response to sinking satellites is still 
beyond the capacity of current AT-body simulations, and so 



the situation remains unclear. We argue, however, that our 
results are qualitatively correct, because the dark matter 
density is the highest in the inner region and so a large frac- 
tion of the orbital energy must be deposited there. Further- 
more, if the mass fraction in cold clouds is as high as ~ 10% 
of the total halo mass, the orbital energy of the clouds be- 
comes comparable to the total binding energy of the dark 
matter in the inner region at a radius r ~ O.lrh- Thus, even 
based on simple energy consideration, the 'heating' of dark 
matter by the cold clouds is expected to be important in the 
inner region. 



3.2 Gas outflow and its effect on the halo profile 

As the dark matter particles gain energy from gcis clouds, the 
dark halo tends to expand to establish a new equilibrium. 
On the other hand, the accumulation of gas clouds in the 
halo centre deepens the gravitational potential, which can 
cause the halo to contract. Which effect dominates depends 
on the details of the system in consideration. However, the 
rapid collapse of cold gas in the halo centre is likely to be 
accompanied by rapid star formation, which may drive a 
large amount of the gas out from the halo centre. Observa- 
tions of nearby starbursts indicate that the mass-loss rate 
in a system is typically 1 to 5 times the corresponding star 
formation rate (e.g. Martin 1999). Such mass loss reduces 
the depth of the potential well and causes the halo to re- 
expand, reducing the net contraction of the halo due to gas 
accumulation. In the extreme case where all the cold gas 
that has sunk into the halo centre due to dynamical friction 
can be expelled from the halo centre as outflow, the change 
of the halo profile due to the accumulation of gas should 
be zero (assuming adiabatic process), and so the only effect 
that can change the halo profile is the energy tranfer from 
the gas clouds due to dynamical friction. Thus, in the sce- 
nario we are considering here, the net effect is always for the 
dark halo to expand. 

If we assume each mass shell expands adiabatically after 
the mass loss, in principle we can solve for the final density 
profile from the hydrostatic equilibrium using phase space 
density profile obtained above. Here, we use a simple den- 
sity model that approximately reproduces the phase-space 
density profile. We assume that the final density profile can 
be described by the form: 



p{r) 



pore 



(rc + r)3 ' 



(27) 



where po is the density at the centre, and Tc is a core radius. 
This profile matches well the NFW profile in the outer re- 
gion. Assuming isotropic velocity dispersion, the hydrostatic 
equilibrium equation is 



rf(pg') _ G'[Mh(< r) + Ma] 
dr ^ 



(28) 



where Mh(< r) is the halo mass within radius r, and Ma is 
the added gas mass (i.e. the total gas mass collapsed into 
the halo centre minus the gas mass ejected by the outflow). 
Together with the adiabatic equation of state p{r)/a^{r) = 
Q(ri), where n is the radius that encloses the mass Mh(< r) 
in the initial density profile, this equation can be integrated 
to give 
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Qin) 
Qh 



p_ 



5/2 r 



p(r) V^(r) dr 



-1 -3/2 



Ph 



(29) 



Here ph and ah are the density and velocity dispersion at 
the virial radius, Qh = ph/e^h, and Vc{r) is the circular ve- 
locity of the halo at radius r. Note that for a given initial 
density profile, ri is equivalent to Mh(< r). With the den- 
sity profile (27), and assuming the value of Qh to be the 
same as that for an NFW profile, we can obtain the phase- 
space density profile Q(ri), or Q[Mh(< r)]. The dashed lines 
in Fig. 4 show the phase-space density profiles obtained by 
assuming different values of c' = rh/rc. For simplicity, we 
have assumed that Ma is negligibly small. As one can see, 
over a large range of mass, the phase space density profiles 
obtained from §3.1 can be matched by the density profile 
(27) with appropriate choices of the value of c = rh/rc. For 
/gas = 0.16, 0.08, 0.04, and 0.02, we obtain c ^ 2.2, 3.5, 5.5 
and 8.0, respectively. The density profile (eq. 27) does not 
match the phase-space density in the inner-most part very 
well, which implies the profile corresponding to the modified 
phase-space density distribution decreases with radius in the 
inner-most region. This is clearly not physical. In reality, the 
inner paxt of the halo must relax to am equilibrium configu- 
ration that depends on the details of the processes involved. 
However, the mismatch involves only a very small fraction 
of the total mass: ^ 2% for an NFW profile with c = 11 and 
even smaller for shallower profiles. We therefore believe that 
the simple model adopted here catches the essence of the 
problem. Fig. 5 shows the rotation curves given by (27) with 
these concentrations, as compared with the original NFW 
profile with c = 4. Thus, depending on the mass of the gas 
that can collapse into the halo centre and ejected by galaxy 
wind, the resulting rotation curve can be flattened substan- 
tially relative to that given by the original NFW profile. 
Within the original scale radius r^, the mass decrease by 
dynamical friction is comparable to the total gas mass in 
clouds. 



4 OBSERVATIONAL CONSEQUENCES 
4.1 The Tully-Fisher relation 

In the slow accretion phase after the formation of the core of 
potential, the virial radius increase roughly as H~^{z) with- 
out affecting much the inner structure. The halo concentra- 
tion will then increase by a factor of about H{zf)/Ho from 
the formation time to the present time. For galaxy halos 
with typical formation redshift ~ 2, this is a factor of 2 to 3. 
Thus, for a pre-processed halo with concentration c' ~ 3.5 at 
formation, the concentration at the present time will be be- 
tween 7 and 10. As shown in Mo & Mao (2000), disk galaxies 
formed in halos with profile (27) and with such concentra- 
tions can match the observed Tully-Fisher relation. Thus 
the flattening of dark matter halos in the fast collapse phase 
by the baryon component can help to reconcile the standard 
ACDM model with the observed Tully-Fisher relation. 



4.2 Dark matter in galaxies 

The concentration of galaxy halos can be studied by looking 
at the fraction of dark matter mass in galaxies. The mass in 




Figure 5. The rotation curves of the NFW profile (solid curve) 
with c = 4 and of the profile (27) with c' = r^/rc = 2.2, 3.5, 
5.5, 8.0 (from bottom up). The rotation velocity is in units of Vh 
while the radius is in units of rh. 



the central parts of galaxies are expected to be more domi- 
nated by daxk matter, if dark halos have higher concentra- 
tions, and vice versa. Since the total mass can be inferred 
from observations of kinematics, while the contribution from 
stars and gas can be inferred by other means, it is possible 
to infer what fraction of the gravitational mass is in dark 
matter. For spiral galaxies, this can be done by modeling 
the rotation curves and the stellar mass distribution in de- 
tail. The results arc not yet conclusive (see e.g. Bosma 1998 
and references therein). The rotation curves for some bright 
spirals are observed to decrease just beyond their optical 
radius, which is best explained if the mass in the central 
region is dominated by the stellar component (e.g. Noorder- 
meer et al. 2003). On the other hand, there are also spiral 
galaxies, whose rotation curves are quite flat around the op- 
tical radius (e.g. Button et al. 2003). The rotation curves 
of many bright spiral galaxies can be fitted by the maximal 
disk model, in which the inner rotation curve is assumed 
to be completely dominated by the stellar components. Al- 
though this does not mean that real spirals all have maximal 
disks, it does suggest that the inner rotation curves of some 
bright galaxies may be dominated by the stellar component. 
Similar constraints may be obtained from the fact that the 
bars in some spiral galaxies are fast rotating. This is not 
expected if the central region of a barred galaxy were dom- 
inated by dark matter, because the bar rotation would be 
slowed down rather rapidly by the dynamical friction of the 
dark matter (e.g. Debattista & Sellwood 2000; Weiner, Sell- 
wood & Williams 2001). All these observations suggest that 
some spirals are dominated by the stellar component in their 
inner regions. 

If the halos of spiral galaxies are as concentrated as 
CDM halos given by A''-body simulations, most spiral galax- 
ies will be dominated by daxk matter, unless the disk is 
unreasonably heavy (e.g. Mo, Mao & White 1998). On the 
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other hand, for a pre-processed halo with an extended core, 
the contribution from dark halo is significantly reduced, 
making it more likely to observe spirals with inner rotation 
curves dominated by stellar mass. 

For elliptical galaxies with moderate luminosities {L ~ 
I/*), recent observations based on the kinematics of plane- 
tary nebulae show that the observed velocity dispersion pro- 
files up to ~ 5 effective radii are consistent with being totally 
due to stellar mass (Romanowsky et al. 2003) , indicating the 
presence of little or no dark matter at these radii. As dis- 
cussed by these authors, this result is difficult to understand 
if these galaxies possess dark halos with properties similar to 
those of the CDM halos. However, if the formation of these 
galaxies is similar to that of the bulges of spiral galaxies, 
their halos may be significantly fiattened by the processes 
envisaged in this paper. In this case, the observational re- 
sult may be easier to understand with the standard ACDM 
model. 



4.3 Survival of substructures in the halo 

After the formation of the core of potential, halo will con- 
tinue to accrete mass in the form of small halos. High- 
resolution Ai'-body simulations show that some of the ac- 
creted halos can survive as subhalos (e.g. Klypin et al. 1999; 
Moore et al. 1999; De Lucia et al. 2003). Intriguingly, the 
predicted number of sub-haloes in CDM exceeds the ob- 
served number of luminous satellite galaxies in a Milky- Way 
type galaxy. One solution for this crisis may be that some 
of the subhalos, especially those of lowest mass, do not form 
many stars and hence remain dark. One of the best ways 
to detect such substructures is through gravitational lens- 
ing. The "anomalous" flux ratios in gravitational lenses are 
thought to be evidence for substructures (e.g. Mao & Schnei- 
der 1998; Kochanek & Dalai 2003 and references therein). 

In an NFW halo, most of the subhalos are distributed 
in the outer part, because subhalos can be tidally destroyed 
by tidal force in the inner part. The situation is different in 
our scenario, where subhalos may exist in the central part 
of a galaxy halo because of the reduced tidal effect in the 
extended core. 

Consider the accretion of a minihalo with circular ve- 
locity ~ 30kms~^ by a galaxy halo with a circular velocity 
of 200kms~^. In the presence of a strong UV background, 
minihalos are not expected to trap gas to form stars (e.g. 
Gnedin 2000), and so their density profiles are not altered 
significantly by the interaction with the baryonic compo- 
nent. If such halos can be approximated as NFW profiles, 
then according to equation (8), they will have a typical con- 
centration c ~ 20. We calculate the tidal radius of such a 
minihalo in a galaxy halo with or without a core. The tidal 
radius, rj , of a minihalo with mass m moving at a radius R 
from the centre of a large halo with mass M is the minimum 
of the following two radii: (1) a radius at which the gravity 
of the small halo is equal to the tidal force of the big halo, 
and (2) a radius defined by the resonances between the force 
the minihalo exerts on the particle and the tidal force by the 
large halo (e.g. Klypin et al. 1999). The first radius can be 
solved from 
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Figure 6. The ratio between the mass within the tidal radius 
and the original virial mass for a small halo (approximated by an 
NFW profile with c = 20) at a given distance from the centre of a 
large halo. The large halo is either assumed to have a NFW profile 
with c = 10, or has the profile given by (27), with c' = ri^/rc = 5, 
10, and 20, respectively. 



where M{R) is the mass within radius R. The second ra- 
dius, assuming primary resonance, is given by r2(r)|sniaii = 
fi(^?)|iarge, wherc n is the angular speed. Fig. 6 shows the 
mass a minihalo can retain within the tidal radius as a func- 
tion of the distance to the centre of the big halo. Results 
are shown for cases where the large halo either has a NFW 
profile with c = 10, or has the profile given by (27) with 
rh/rc ~ 5, 10, and 20, respectively. If the halo has a core with 
c' = rh/rc as small as 10, as is required to match the Tully- 
Fisher relation, a minihalo can retain more than 20% of its 
original mass, in strong contrast with the case where the 
host halo has the original NFW profile. This clearly demon- 
strates that the survival of subhaloes in the central part of 
galaxies is likely to be sensitive to the inner mass profiles. 
Note, however, that in Fig. 5, we have neglected the effect 
of baryons which may steepen the central density profiles 
and hence destroy some subhaloes. It is currently unclear 
what might be the eventual survival fraction of subhaloes in 
a realistic galaxy. 



5 PREHEATING AND GALAXY FORMATION 

5.1 Preheating of a protogalaxy region 

In order for the starburst in the core of potential to eject a 
large amount of gas from the dark halo potential well, suffi- 
cient number of stars must form so that the kinetic energy 
from supernova explosions from massive stars can heat the 
gas to a temperature comparable to the escape temperature. 
For an NFW profile, the escape velocity can be calculated 
from equation (7). We obtain 
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(31) 



where /3 ^ 3.1 and 3.7 for c = 4 and 10, respectively. We 
write the star formation rate as 



Ttff 



(32) 



where Mgas is the total mass of cold gas that has assembled 
in the halo centre in the fast accretion phase, iff is the free- 
fall time of the gas, and t is a constant. The total mass of 
stars that can form can be written as 



(33) 



where is the star formation time, and e* = tt,/{Ttff) is 
the star formation efficiency. The total kinetic energy from 
Type II supernova explosions can be written as 



(34) 



where v is the number of supernovae per unit mass, and esn 
is the kinetic energy per supernova. We adopt esn = 10®^ erg 
and V = (125 M©)"^. The value of v quoted here assumes a 
Salpeter initial mass function (with a lower cutoff of O.IM© 
and an upper cutoff of lOOM©). If a fraction of eo of this 
energy is to heat the gas, the specific thermal energy of the 
gas will be 



(630 km s ^)^eoe*- 



(35) 



Requiring this specific energy to be equal to the escape tem- 
perature, we obtain the star formation efficiency. 

Note that this relation is meaningful only when e* < 1, and 
so it only applies for halos with circular velocities smaller 
than e(y^/3^^630kms~^. For larger halos, all cooled gets can 
form stars and so 1. This scaling relation between star 
formation efficiency and halo circular velocity is similar to 
that obtained by Dekel & Silk (1986) based on a calculation 
of supernova remnants in a uniform medium, and has been 
extensively used in semi-analytic models of galaxy forma- 
tion (e.g. White & Frcnk 1991). The value of eo obtained 
by Dekel and Silk is about 0.02. Since e* must be smaller 
than 1, this low value of eo implies that wind can only es- 
cape in small halos, with Vh ^ 30kms~^ (assuming /3 = 3). 
On the other hand, observations of supernova-driven wind 
in starburst galaxies require much higher efficiency. The ob- 
servations of superwind in local starburst galaxies suggest 
that most of the supernova energy is contained in the su- 
perwind (Heckman et al. 2000). Such high efficiency may 
be more relevant to our discussion, as we arc mainly con- 
cerned with starbursts. In this case, wind can escape from 
halos with circular velocities as high as 200 km s"'^. In prin- 
ciple, the effective value of eo can be larger than 1, if active 
galaxy nuclei associated with starbursts also contribute to 
the heating. 

The gas that has been heated up to the escaping tem- 
perature will escape the core of potential on a time scale 
teac ~ r-cp/Ksc ^ l/[W0H{z{)], whcrc rep is the radius of 
the core of potential, and H(zf) is the Hubble constant at 
the redshift zt when the core of potential was formed. This 
time scale is shorter than the time scale for mass accretion 
in the slow accretion phase, which is typically the age of the 



universe. Thus, the out-going wind will interact with the 
gas outside the core of potential before the gas is accreted 
by the potential well. The average overdensity within a re- 
gion that will eventually form a halo at the present time 
is 5{zf) = 5c/D{zf), where 5c ^ 1.68 is the critical linear 
overdensity for spherical collapse, and D{zf) is the linear 
grow factor at Zf. For a flat universe with Qm.o = 0.3 we 
have S{zi) fs 0.43, 0.53, 0.71, 1.03 for zi = 4, 3, 2 and 1, re- 
spectively. Thus, for halos with an early formation time, the 
overdensity around the core of potential is moderate at Zf. 
Bcised on the extended Press-Schechter formalism (Bond et 
al. 1991), one can show that, in the standard ACDM model, 
more than half of the mass of a present galaxy-sized halo 
was in progenitors with circular velocities ^ 50kms~^ at 
redshift z = 2 - 3. Such progenitors can retain only a small 
fraction of their gas either because their potentials are too 
shallow to trap gas heated by photoionization, or because 
supernova explosions in themselves or from nearby galax- 
ies can effectively expel the gas (e.g. Scannapieco, Ferrara 
& Broadhurst 2000). Thus, at the formation redshift of a 
present-day galaxy halo, most of the mass that is not in the 
core of potential is in halos that cannot trap much gas. Thus, 
the wind generated in the core of potential by starburst will 
be propagating in a relatively uniform medium. 

Because the gas density in the post shock medium is in 
general quite low and so radiative cooling of the shocked gas 
is negligible, we may use the wind solution of Weaver et al. 
(1977) (which assumes a constant energy injection rate and 
neglects radiative cooling and galaxy potential) to get some 
idea about how much gas can be affected by the wind. In 
this case, the radius of the shock front changes with time as 



/ 125 / EAf 



1/5 



(37) 



where E is the total supernova energy, p is the gas density 
of the medium, and At is the propagation time. Inserting 
E = (l/2)M„ind/3^V'h into the above expression, and using 
the relations between Mh, rh and Vh, we obtain 



Ra ~ lOr-h 





f-windP 


1/5 
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(38) 



where tf m 1/H{z) is the age of the universe when the 
starburst is generated, which is set to equal the formation 
time of the core of potential, and /wind is the gas mass 
fraction in the wind. This radius corresponds to a mass 
M ~ (1000/Ah)(Af/t/)®/^Mh, where Mh is the mass of the 
core of potential, i.e. the mass of the halo at the formation 
time. Thus, within a Hubble time, the wind can affect a 
large volume around the halo. For a halo that formed at a 
redshift Zf with circular velocity Vh and concentration c, the 
total mass at the present time is about H{zi)/Ho times the 
progenitor mass. Thus, for a halo formed at 2 ~ 3, all the 
gas in the halo proper can be affected by the wind generated 
in the core of potential, if the star formation efficiency is as 
high as discussed above. 

Unfortunately, the details about the interaction be- 
tween the outgoing winds and the intergalactic medium are 
very complicated, and wc arc not able to make a detailed 
discussion here. If we assume that half of the energy of the 
wind is thermalized in the protogalaxy region, then the spe- 
cific thermal energy gained by the gas is 
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£• = [Ho/H(zf)](/3\/h)V4. 

The corresponding specific entropy of the gas is 

[g0/Jj-(Zf)]/wind(/3H)' 



(39) 



5 = 



T 



4[n(0f)(l+5/)]2/3 



8.5 X lO^keVcm^ 



lHo/H{zi) 



(l+2f)2(l + 5^)2/3 



200 km s- 



(40) 



where n(zf) is the mean gas density of the universe at Zf, 
and 5/ is the overdensity of the protogalaxy region. Thus, 
in terms of the value of S, heating is more effective for ha- 
los with higher Vh and lower zj. Taking Z{ = 3, 5f — 1, 
and /3 = 3, we have S ~ 75(14/200 kms-i)2keVcm2. 
This should be compared with the characteristic specific 
entropy of the gas in a virialized halo: Si, = T^jruJ^ ~ 
70keVcm^(Vh/200kms-i)V(l + «)^- Note that both S and 
Sv scale with Vh in a similar way, and S ^ Sv at all z ^ 0. 
Thus the heating from the core of potential can have an im- 
portant impact on the subsequent collapse of the gas into 
galaxy halos. 



5.2 Formation of disk galaxies 

In the scenario we are considering here, it is quite natural 
to associate the formation of galaxy bulges with the fast 
collapse phase. Depending on whether a significant disk can 
grow in the subsequent slow accretion phase, a spiral galaxy 
with some disk/bulge ratio will form. Since the efficiency 
of gas cooling is reduced in a preheated gas, the amount of 
gas that can eventually form a disk can be much lower than 
that implied by the universal baryon fraction (see Mo & Mao 
2002 for details). Also, the bulge/disk ratio is expected to 
depend on environment , because galaxy halos in high density 
regions may not be able to accrete much gas before they 
merge into large systems where radiative cooling is no longer 
effective. As pointed out in Mo & Mao (2002), this may 
be responsible for the observed morphology segregations of 
galaxies. In what follows, we consider several other factors 
that may have played important roles in determining the 
morphological types of galaxies. 

As shown by equation (36), the star formation efficiency 
in the fast collapse phase is likely to be more efficient in halos 
with higher circular velocity. This effect alone would mean 
that bulges are more important in systems with larger Vh (or 
mass). This is consistent with the observational trend that 
later typo galaxies generally have lower luminosities (e.g. 
Roberts & Haynes 1994). If the halo is massive enough, this 
efficient star formation may heat all the gas around the halo 
to such a high temperature that further cooling of the gas is 
prohibited, making it impossible for the formation of a disk 
component. This might bo the reason why there appears to 
be an upper limit on the masses of disk galaxies (e.g. Peebles 
& Silk 1990). 

In addition to this mass sequence, halo formation his- 
tory may also play a role in determining the bulge-to-disk 
ratio of the galaxy that forms in the halo. Although halos 
of a given mass have a typical formation time, the disper- 
sion among different halos is large for low-mass halos, and 
so there are significant number of present-day galaxy halos 



that can have their cores of potential formed only recently 
(e.g. Zhao et al., 2003b). Such halos in general have low 
concentration, as discussed above, but higher specific angu- 
lar momentum (e.g. Mailer, Dekel & Somerville 2002). Gas 
accretion history by such halos is expected to be different 
from that for a halo with early formation. First, because of 
the late formation, the gas associated with the region that 
eventually collapses to form the core of potential may have 
already been heated up by star formation in small progen- 
itors before the formation of the core of potential. Second, 
since the characteristic density is lower at lower redshift, 
cooling is less effective at the formation time for a later for- 
mation. Thus, the formation of the core of potential of such 
halos is not expected to be associated with rapid formation 
of cold clouds and with strong starburst. Such systems are 
therefore expected to have small bulges and may be iden- 
tified as late-type galaxies. The halos of such galaxies are 
therefore expected to have lower concentration, and the in- 
ner profiles of their halos are expected to be less affected 
by the bulge formation. Because of the late formation, these 
systems are expected to show weaker clustering than nor- 
mal galaxies (e.g. Mo & White 1996). The latter prediction 
is consistent with the observations that late-type galaxies, 
such as low surface-brightness galaxies, have weaker correla- 
tion in space than average spiral galaxies (e.g. Mo, McGaugh 
& Bothun 1994). 

An interesting question here is whether there are truly 
bulgeless galaxies. The halos of such galaxies should not have 
been pre-processed by the bulge formation and so should 
preserve their original profile. Any bulgeless galaxies with 
halos that contain extended cores should therefore be con- 
sidered as serious evidence against the CDM cosmogony. As 
mentioned earlier, the rotation curves of many low surface- 
brightness galaxies are better fit with halo profiles with 
cores. Whether or not all these galaxies have a significant 
bulge component is still a question of debate. Giant bright 
low surface-brightness (LSB) disk galaxies such as MALIN- 
1 are known to have normal bulges (e.g. Bothun, Impey & 
McGaugh 1997) and so their halos are expected to be signifi- 
cantly pre-processed. Early observations in blue bands show 
that many low- luminosity LSBs are bulgeless (see Bothun 
et al. 1997), but recent observations of such galaxies in near 
infrared bands suggest that many of them have detectable 
bulges in old stars and show episodic star formation (e.g. 
Galaz et al. 2003). Since these galaxies are associated with 
shallow potential wells, the implied amount of early star for- 
mation may be significant in driving large amount of gas out 
from halo centre, thereby altering the inner profiles of their 
halos. Clearly, detailed modeling is required in order to show 
if the scenario we are proposing here is consistent with this 
population of galaxies. 



6 DISCUSSION 

We have outlined a scenario of galaxy formation in which 
the host halos of present-day galaxies may be significantly 
flattened in the inner region by interaction with the gas com- 
ponent in an early phase of rapid collapse accompanied by 
starburst. Our scenario is based on recent simulation results 
that such a phase is expected in the current ACDM model of 
structure formation. We have shown that this scenario can 
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help to alleviate several vexing problems in current theory 
of galaxy formation. 

Much further work is required to make quantitative pre- 
dictions with the scenario we have proposed. For example, 
the treatment of dynamical friction is performed in an over- 
simplified manner. In order to model this process quantita- 
tively, high-resolution numerical simulations are required. In 
addition, the assumption of adiabatic contraction should be 
more carefully checked, particularly for cored halo profiles 
(such as that described by eq. 27). Ideally, we should use 
cosmological simulations that fully account for all the rel- 
evant processes, such as gas cooling, cloud formation, star 
formation, outflow and dynamical friction. Unfortunately, 
such simulations arc not yet feasible at the present time. 
A more realistic approach is to tackle the problem step by 
step. For example, the dynamical friction may be studied 
using high-resolution iV-body simulations where cold clouds 
are represented by simple A^-body systems; the interaction 
between the starburst-drivcn wind and the protogalaxy gas 
may be studied using coutrolod simulations of individual sys- 
tems. We intend to return to some of these issues in future 
papers. 
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